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We report on further progress in our programme of understanding confinement in 3d and 4d SU(2) gauge theory 
in terms of Z(2) monopoles. A sufficient condition for confinement was previously translated into Z(2) monopole 
correlation inequalities in a related SO (3) gauge theory. We shall discuss the physical picture underlying this 
scenario and present some Monte Carlo evidence concerning the monopole correlation inequalities. 
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It is an old idea that confinement in Z(N) gauge 
theories can be understood in terms of Z(N) vor- 
tices (also called fluxons) linking with the Wil- 
son loops (see e.g. [0). In three dimension^ a 
Z(N) vortex is a closed stack of plaquettes (form- 
ing a loop) with a nontrivial Z(N) element on 
them. These are excitations very much analogous 
to Peierls contours in the 2d Ising model. 

This vortex condensation mechanism was gen- 
eralised to SU(N) gauge theories by Mack and 
Petkova ||]. Here the magnetic vortices also be- 
long to Z(N), the centre of SU(N) but there is 
a major difference as compared to the Abelian 
case. While the Abelian magnetic fluxon neces- 
sarily has a thickness of one lattice spacing, in 
SU(N) lattice gauge theory, due to the continu- 
ous nature of the group, it can be spread out over 
several lattice spacings. In fact, it is exactly the 
spreading of the vortices that makes it possible for 
them to survive at weak coupling. It was empha- 
sized by several authors that the sufficiently fast 
(exponential) spreading of the magnetic flux is 
an essential feature of any confining SU(N) gauge 
theory ||^,^. The general picture emerging from 
these investigations was that it is probably the 
thick, spread out vortices interlocking with the 
Wilson loop that are responsible for confinement 
at weak coupling. Unfortunately it is extremely 
hard to trace the dynamics of these fluxons at 
different length scales and although they provide 
a nice intuitive picture, not much progress was 
made in obtaining a quantitative understanding 



^Throughout this paper we shall consider gauge theories in 
three dimensions but our arguments can be easily adopted 
to the physically more relevant 4d case. 



of the mechanism. 

For the past few years there has been a 
programme to achieve this goal and rigorously 
demonstrate that confinement persists down to 
arbitrarily small couplings in SU(2) gauge the- 
ory In the following we shall summarise the 
main ideas underlying this programme. Here we 
can present only a rough intuitive picture and for 
more technical details the reader is referred to jj] 
and especially |^ . The latter is the most detailed 
presentation so far, although in the context of the 
analogous 2d SU(2) principal chiral model. 

In the original formulation of Mack and 
Petkova, thick vortices were detected by intro- 
ducing a vortex container, special boundary con- 
ditions along a thick torus linking with the Wil- 
son loop, to trap the vortex inside the torus. In- 
stead of this we essentially use the whole lattice 
as a vortex container. The Wilson loop confine- 
ment criterion can be translated into an exponen- 
tially fast spreading of a vortex winding around 
the cubic lattice, as a function of the finite lattice 
size Q. In order to detect thick vortices winding 
around the lattice we use an equivalent represen- 
tation of the SU(2) gauge theory in terms of a 
dual Ising model (Z(2) gauge theory in the 4d 
case) interacting with an SU(2)/Z(2) gauge the- 
ory. Since the action of the SU(2)/Z(2) system is 
insensitive to the Z(2) centre, a thin Z(2) vortex 
(i.e. a stack of plaquettes) winding around the 
lattice does not cost an energy proportional to 
its length, however due to the periodic boundary 
conditions it signals the presence of a thick vortex 
winding around the lattice. In this way in order 
to trace thick vortices we have to detect only thin 
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ones. 

Contrary to Z(N) gauge theories, SU(N) gauge 
theories admit Z(N) monopoles, elementary lat- 
tice cubes that have a net Z(N) flux flowing out of 
them. (In the Z(N) case these are absent due to 
the Bianchi identity.) These monopoles are gauge 
invariant objects, not like the U(l) monopoles 
that appear in this model only after a partial 
gauge fixing. For a Monte Carlo study involving 
both types of monopoles, see 

The presence of Z(2) monopoles means that in 
our SU(2) model Z(2) vortices can either be closed 
or bounded by a pair of monopoles (a monopole 
loop in 4d). A vortex going all the way around 
the lattice is globally quite similar to the one 
which is bounded by a nearby pair of monopoles; 
the difference between them is only a small local 
excitation. In the Z(2)xSU(2)/Z(2) formulation 
the monopoles of the SU(2)/Z(2) gauge theory 
turn out to couple to the dual Ising model and 
provide an external magnetic field for the Ising 
spins. If they can produce a nonzero effective 
magnetic field for the dual Ising model than the 
Ising spins are ordered which in turn means dis- 
order i.e. confinement in the original gauge the- 
ory. The presence of a nonzero effective external 
field for the dual Ising model can be rigorously 
established provided that two types of monopole 
correlation inequalities are satisfied. 

At this point we want to emphasize that the 
presence of Z(2) monopoles is not necessary for 
confinement. It is of course not small monopoles 
but large spread out vortices that are really re- 
sponsible for confinement. However, as we have 
already seen, some of the large vortices naturally 
occur "tagged" with a pair of monopoles, i.e. they 
are not closed but have a small gap instead. It is 
exactly these tagged vortices that couple to the 
dual Z(2) system and make it possible to estab- 
lish confinement in terms of the Z(2) variables. 
In other words, if we were to constrain out Z(2) 
monopoles, confinement would not be lost but we 
would not be able to see it in terms of the Z(2) 
variables. 

Now returning to the monopole inequalities, 
the first of these is a factorisation inequality per- 
taining to monopole pairs (loops in 4d), this 
was discussed to some length in [Q. The sec- 



ond inequality that we would need to prove is 
that the probability of having a nearby monopole 
pair with their flux winding around the lattice, 
is greater than some finite non-zero number. In 
other words, the free energy of the above config- 
uration is a bounded function of the lattice size. 
A simple semiclassical estimate supports this as- 
sumption in 3 and 4 dimensions (but not above 4). 
In the remainder of the present paper we briefly 
present the results of a Monte Carlo measurement 
of this quantity as a function of the lattice size. 

We could do the measurement only in the 3d 
case since in 4d due to the long lifetime of the 
large vortices we could not obtain enough statis- 
tics. We measured the probability of the config- 
uration with two fixed adjecent monopoles, their 
flux winding around the lattice; normalised by the 
probability of having no monopoles at all and no 
vortex going around the lattice in the given direc- 
tion. The measurement was done by generating a 
series of configurations using a local heat bath al- 
gorithm with the Villain form of the SU(2)/Z(2) 
action 0, 

5[C/,a] =/3^(7ptrC/p, (1) 
p 

where the cr^'s are Z(2) valued plaquette variables 
that are summed over in the partition function 
to ensure the desired Z(2) invariance. We had 
to ensure that the P we chose was already in the 
weak coupling region since we were interested in 
the weak coupling behaviour of the quantity mea- 
sured. Unfortunately in 3d the crossover between 
weak and strong coupling is not so sharp as it 
is in 4d and in 3d the specific heat peak at the 
crossover is missing. Therefore we used the expo- 
nential falling of the monopole density as a cri- 
terion for weak coupling. We chose (3 — 11.5, as 
can be seen in Figure this is already well in the 
exponential regime of the monopole density. 

Since (at a given /3) on larger lattices there 
are typically more monopoles, the probability 
of the confugurations with exactly one an zero 
monopole pairs decreased very rapidly with in- 
creasing lattice size. This meant that the quan- 
tity we measured was given as a ratio of two num- 
bers, both becoming very small on larger lattices, 
however their ratio was expected to be fairly sta- 
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Figure 1. The Z(2) monopole density as a func- 
tion of /3 on a 8'^ lattice. 



ble. This made the signal less accurate on larger 
lattices and eventually prevented us from going 
beyond a lattice size of 10'^. Even at this point we 
typically needed several thousand configurations 
to get a signal at all. Our results are summarised 
in Figure ^ which shows the measured probabil- 
ity as a function of the lattice size. As we expect, 
the probability of the two- monopole configuration 
with the long fluxon does not decrease with the 
lattice size. This is consistent with the semiclas- 
sical estimate and makes it quite improbable that 
it can go to zero as L — > cxd. 

If confinement is present we expect an exponen- 
tial spreading of the flux, i.e. its free energy to 
approach zero exponentially on sufficiently large 
lattices. In our framework we would need to prove 
only the substantially weaker condition that the 
flux free energy does not diverge with the lattice 
size. It is remarkable that even this weaker prop- 
erty is absolutely nontrivial to prove rigorously 
and we had to resort to Monte Carlo. 

Finally we would like to point out that each 
step in the prgramme outlined here has been 
given precise formulation and to complete it we 
would need to prove the correlation inequalities 
discussed above. 
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Figure 2. The probability of two adjecent 
monopoles with their connecting fluxon going 
around the lattice divided by the probability of 
having no monopoles at all and no fluxon going 
around the lattice in the given direction. 
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